Background reionization history from omniscopes 
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The measurements of the 21 cm brightness temperature fluctuations from the neutral hydrogen at 
the epoch of reionization (EoR) should inaugurate the next generation of cosmological observables. 
In this respect, many works have concentrated on the disambiguation of the cosmological signals 
from the dominant reionization foregrounds. However, even after perfect foregrounds removal, our 
ignorance on the background reionization history can significantly affect the cosmological parameter 
estimation. In particular, the interdependence between the hydrogen ionized fraction, the baryon 
density and the optical depth to the redshift of observation induce non-trivial degeneracies between 
the cosmological parameters that have not been considered so far. Using a simple, but consistent, 
reionization model, we revisit their expected constraints for a futuristic giant 21 cm omniscope by 
using for the first time Monte-Carlo-Markov-Chains (MCMC) methods on multi-redshift full sky 
simulated data. Our results agree well with the usual Fisher matrix analysis on the three-dimensional 
flat sky power spectrum but only when the above-mentioned degeneracies are kept under control. 
In the opposite situation, Fisher results can be inaccurate. We show that these conditions can 
be fulfilled by combining Cosmic Microwave Background measurements with multiple observation 
redshifts probing the beginning of EoR. This allows a precise reconstruction of the total optical 
depth, reionization duration and maximal spin temperature. Finally, we discuss the robustness of 
these results in presence of unresolved ionizing sources. Although most of the standard cosmological 
parameters remain weakly affected, we find a significant degradation of the background reionization 
parameter estimation in presence of nuisance ionizing sources. 

PACS numbers: 98.80.Cq, 98.70.Vc 



I. INTRODUCTION 

Among the next generation of cosmological probes, in- 
terferometric radio telescopes observing the redshifted 
21 cm line associated with the hyperfine transitions of 
neutral hydrogen atoms have attracted a lot of attention 
(see Refs. pj-j4|] for reviews). These telescopes are revo- 
lutionary in their design as in the proposed Fast Fourier 
Transform Telescope (FFTT) which is conceptually an 
all digital antennas array imaging the whole visible sky 
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at once 'S!\. Images would then be reconstructed by a 
two-dimensional fast Fourier transform over the N an- 
tennas signal, i.e. in JVlniV operations. As argued in 
Ref. @, compared to the required A^ 2 pairing in tradi- 
tional interferometers, the gain could be used to scale up 
the telescope size and sensitivity thereby allowing mea- 
surements of cosmological signals. The first generation 
of such telescopes is under deployment [714131] . Although 
not designed for cosmological purposes, they aim at de- 
tectinga cosmological signal from the epoch of reioniza- 
tion [1J-[20|. It is still a matter of active research to 
know if the signal coming from the astrophysical fore- 
grounds can be properly separated from the cosmologi- 
cal one [21-25]. One should indeed keep in mind that 
the former is actually a few order of magnitude stronger 
than the latter and this has triggered interests into the 
cross-correlation of the 21 cm signal with much cleaner 
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FIG. 1: Background evolution of the neutral hydrogen gas 
and spin temperatures (top) as a function of the redshift dur- 
ing the dark ages. The spin temperature T s being defined 
as m/no = 3 exp[— FtW2\ / (kBT a )\ where no and m denotes 
the number density of atoms respectively in the singlet and 
triplet hyperfine states. The bottom panel shows the resulting 
brightness temperature. 



data such as the Cosmic Microwave Background (CMB) 
or galaxy surveys [2rjl430l |. 

The physical origin of the cosmological 21 cm radiation 
lies in the differential cooling induced by the expansion 
of the universe on relativistic and non-relativistic gases. 
After recombination, one would naively expects the tem- 
perature of neutral hydrogen T gas to decrease in I /a 2 , 
where a is the Friedmann-Lemaitre-Robertson- Walker 
(FLRW) scale factor, which is faster than the radiation 
temperature T ra d scaling as I /a. As a result, the neutral 
hydrogen gas becomes cool enough to be able to absorb 
CMB photons by a spin flip hyperfine transition at a 
wavelength of 21 cm. Tuning a radio-telescope at the 
corresponding rcdshiftcd frequency allows, in principle, 
to probe the density fluctuations of neutral hydrogen over 
the sky and at an y re dshift, hence the so-called cosmolog- 
ical tomography [3l|, [32( ■ In reality, the situation is a bit 
more complex and one has to take into account the evo- 
lution of the Boltzmann distribution of neutral hydrogen 
due to the collisions with the residual electrons, protons 
and absorption versus stimulated emission of CMB pho- 
tons [33 3G] . The background evolution can nevertheless 
be numerically derived and we have plotted in Fig. [T] the 
resulting brightness temperature evolution. 

For an assumed generic set of cosmological parameters, 
taken from the Wilkinson Microwave Anisotropics Probe 
(WMAP) seven year data [13] , the absorption is maximal 
around a redshift of z ~ 10 2 . Notice that at lower red- 
shifts, the spin temperature is driven again towards the 
photon temperature and the signal vanishes till the epoch 
of reionization (EoR). The evolution of neutral hydrogen 
density fluctuations can similarly been predicted during 



FIG. 2: Reionization driven spin temperature (top) and 
brightness temperature (bottom) with respect to the observa- 
tion redshift. The reionization model is described in Sec. Ill Bl 
(r = 0.1 
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the dark ages by using the theory of cosmological pertur- 
bations 38- 40]. Provided the length scales of interest re- 
main in the linear regime, the theoretical predictions for 
the 21 cm power spectra are neat (4ll444| . Being three- 
dimensional in nature, the information content is huge 
and have been used to forecast constraints on various 
cosmological models such as non-Gaussianities 
cosmic strings l4^M5li , dark matter signatures 
modified gravity |57| or inflation [58146 lj . 

As already mentioned, in addition to have a small am- 
plitude compared to foregrounds, the dark ages signal 
is redshifted to wavelengths of typically twenty meters 
and it is not obvious that it may actually be used for 
cosmology in a foreseeable future. One very futuristic 
approach would be to build those radio-telescopes on the 
Moon [62] ■ A more reasonable approach is the possibility 
to probe the hydrogen density fluctuations with the 21 
cm line at the epoch of reionization [63M65| . This time, 
the exciting photons are coming from reionization sources 
(stars or quasars) instead of the CMB and are expected 
to produce an efficient spin states population inversion. 
This results into T s 3> T ra( j and generates an emission 
line at a much smaller redshifted wavelength. We have 
plotted the expected evolution of the brightness temper- 
ature during the epoch of reionization in Fig. [51 for a 
simple reionization model that we describe in Sec. Ill Bl 

The resulting signal is around 25 mK in amplitude, 
which is of the same order than those coming from the 
dark ages, but at a higher frequency thereby rendering 
its detection more likely. However, the theoretical pre- 
dictions are now far less neat due to the additional as- 
trophysical uncertainties associated with the way reion- 
ization proceeds [66|, H3]- However since 21cm data 
are three-dimensional by nature, it has been shown in 
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Rcfs. [68, 69] that the redshift evolution could be used to 
efficiently eliminate the expected foregrounds while still 
keeping some of the so-called longitudinal modes for cos- 
mology. 

Forecasting the 21 cm constraints for cosmology is usu- 
ally made from Fisher matrix analysis which merely as- 
sumes the likelihood to be Gaussian around the best 
fit. Provided the parameters are well constrained, the 
method is fast and accurate [lO]. However, as noted 
in Ref. [Hj], some model parameters linked to reioniza- 
tion, such as the optical depth or the hydrogen ionized 
fraction, can be completely degenerated from the 21 cm 
point of view and factorize out of the Fisher matrix. In 
a realistic situation, this should not be the case as all 
reionization parameters are uniquely determined by the 
background reionization history. One may therefore won- 
der how their correlations to the standard cosmological 
parameters affect the forecasts. 

In Refs. [20|, [7lJ, the problem of reconstructing the 
background reionization history is specifically addressed 
in the context of future 21 cm experiments that would 
be devoted to the global signal, i.e. the homogeneous 
mode. Here, we would like to discuss this issue for the 21 
cm FFTT-likc experiments which arc poorly, if not at all, 
sensitive to the constant mode. In that situation, the cos- 
mological brightness fluctuations ST^(x, z) have simulta- 
neously a dependence in both the background evolution 
and spectral shape. This question is peculiar to 21 cm ob- 
servables because the homogeneous mode T\ } (z) depends 
on various cosmological parameters, and when concerned 
with EoR, on the reionization (see Sec lII A|) . Let us notice 
that such a situation does not occur for CMB not only 
because T cm b is measured, but because it is uniquely de- 
termined by Q-y. 

In this context, we would like to quantify how cosmo- 
logical forecasts are affected by the background reion- 
ization history, and if it is possible for omniscopes to 
constrain the reionization parameters. For this purpose, 
we consider a giant FFTT-like ground based instrument 
and assume that most of the foregrounds can indeed be 
eliminated using redshift evolution. As a result, we will 
be keeping only a few redshift slices for our forecasts, 
eventually marginalizing over some nuisance ionization 
spectra modelled as in Ref. [44j. We consider a sim- 



ple reionization model for which the time evolution of 
the spin temperature, hydrogen ionized fraction and op- 
tical depth are given as a function of some extra pa- 
rameters. In this manner, it is imposed that they are 
intrinsically correlated by reionization physics. As al- 
ready mentioned, in order to go beyond Fisher matrix, 
we use MCMC methods over a modified version of the 
CAMB 21 cm code [H, |H, [72[ which includes the epoch of 
reionization as well as some nuisance ionization spectra. 
This allows us to use and compare both Fisher matrix 
analysis and full sky simulated data coupled to MCMC 
methods [73|. 

The paper is organized as follow. In Sec. [ill we present 
the reionization and telescope models used in the analy- 



sis while the forecasts methodology and results are pre- 
sented in Sec. IIIII Starting from one redshift slice only, 
we show that the reionization induced degeneracies dam- 
age, and sometime prevent, the determination of some 
cosmological parameters, such as the baryon density. In 
these situations, we also find that a Fisher analysis can 
be inaccurate for some strongly correlated parameters. 
Combining a few redshifts, which sample the beginning 
of reionization, improves the situation but does not com- 
pete with a Planck-like CMB experiment. However, we 
show in Sec. lIIIC~3l that combining 21 cm and CMB lifts 
almost all those degeneracies thereby improving the over- 
all accuracy by an order of magnitude. In that situation 
our MCMC results match with the usual Fisher matrix 
analysis while the underlying reionization model can be 
completely determined. Finally, we show that these re- 
sults may be mitigated by the presence of toy nuisance 
reionisation spectra. 



II. REIONIZATION AND TELESCOPE MODELS 
A. 21 cm signal from the EoR 

1. Background 

Assuming a small emission line profile, the 21 cm 
brightness temperature durin g th e EoR and at a given 
observed energy e is given by [42J 



T b (e) = (1-e- 



-*g -*rad 



l + z 

where r e is the optical depth to 21 cm 
3fic 3 A w n m (rj) 



l&k B vl 1 T s {rf)H(rf} 



(1) 



(2) 



All quantities are evaluated at the conformal time r\ = r) e 
at which e = a(j} e )E2i = 2ira(r] e )hv2i, ^21 being the 
rest frame frequency of the 21 cm spin flip atomic tran- 
sition and a is the scale factor. In the following, we 
use natural units fc B = K = c = 1. In Eq. (|2]). 
A w ~ 2.869 x 10~ 15 s^ 1 is the Einstein's coefficient of 
spontaneous emission, n HI (r]) the density of neutral hy- 
drogen atoms and H (77) the Hubble parameter. Introduc- 
ing the neutral fraction x H = n HI /(n HI + Tii), with m the 
number density of ionized hydrogen, one can express the 
density of neutral hydrogen in terms of the cosmological 
parameters today as 

3M 2 (1 



n H1 (z) = x H (z)- 



Y) 



(3) 



S7rG N m H 

where Y ~ 0.24 is the helium mass fraction 1 and m H 
the mass of the hydrogen atom. During reionization, one 



1 Notice that the Helium fraction affects significantly the bright- 
ness temperature for a given Ht> . 
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assumes that T s is driven to high values due to ionizing 
sources [67j such that one can expand Eq. ([TJ in 



T, 



1 



-T ( 



O 



(4) 



where 



T, 



T f T R = 



9A 10 Mg,n b gg(i-y) 

l&v^mn a 3 (7j e )H(r] e ) 



X H (Ve)- (5) 



In the limit T s > T £ , T b ~ aT e . The brightness temper- 
ature is thus positive (emission line), does not depend 
anymore on the spin temperature and is directly propor- 
tional to x H . Up to the metric fluctuations, this is also 
the case for the linear perturbations (5T b which directly 
probe the neutral hydrogen fluctuations. However, the 
approximation T s 3> T e cannot hold at the beginning of 
reionization as T s should continuously rise from the dark 
ages value to its maximum value deep in the reioniza- 
tion era. As a result, it is far from evident that there 
exists a redshift range for which this approximation is 
fully consistent. Actually, even for a spin temperature 
T s ~ 10 3 K, one can check that differences at the percent 
level already arise between Eq. (JTJ and the leading term 
of Eq. (H]). For these reasons, in this paper, we keep the 
complete dependence in T s and T ra d and assume a sim- 
ple, but consistent, model of reionization to parametrize 
the spin temperature evolution (see Fig.[2]and Sec. Ill B[ ). 



2. Linear perturbations 

The brightness temperature fluctuations can be ob- 
tained by solving the Boltzmann equations driving the 
evolution of the 21 cm photons in a perturbed FLRW 
space-time. This has already been done by Lewis and 
Challinor during the dark ages in Ref. [42j ■ We have fol- 
lowed their approach while adding the contributions of 
the hydrogen ionized fraction perturbations coming from 
the reionization sources. For this pur pos e, we have mod- 
ified the publicly available CAMB code [72|, uM an< l numer- 
ically solved the full Boltzmann equations incorporating 
all the linear order effects discussed in Ref. [42j], i.e. ne- 
glecting only anisotropics higher than dipole photon ones 
and any broadening of the emission line profile. The mul- 
tipoles components of the 21 cm distribution function ac- 
tually used in the next sections have been written down in 
Appendix [A] In order to illustrate the physical processes 
at work, it is nevertheless convenient to approximate the 
perturbed Boltzmann distribution in the small scale limit 
and at leading order in t £ . In the line of sight direction 
n, at the measured energy e, one gets 



<5T b (x,n,e)~T b e- T <= 



A 1 ~ dV 

A s n ■ — — 

aH or) 



where the A s is monopole source 

-^rad 



A s = A HI 



(A Ta - A Trad ), 



(6) 



(7) 



and all A x = Sx/x stand for relative perturbations. The 
second term of Eq. © encodes the redshift distortions 
due to the perturbations v in the (baryonic) gas rel- 
ative velocity to the observer. As shown in Ref. [42j], 
Thompson scattering suppresses the 21 cm brightness on 
all scales and is responsible of the exponential term in 
Eq. ((6]), T c (rj e ) being the Thompson optical depth to the 
redshift of observation (not to be confused with r e ). Dur- 
ing the dark ages, x H ~ 1 and the perturbations A m in 
the neutral hydrogen directly trace the baryonic pertur- 
bations Ab = <5n b /n b . At the EoR, the ionizing sources 
on their own are expected to induce new perturbations 
in the neutral fraction such that A m = A^ + At, (at 
fixed Helium fraction), or, in terms of the hydrogen ion- 
ized fraction x t (not to be confused with the total ionized 
fraction x e ) 



A H I 



(8) 



where A Xi is the relative ionized fraction perturbation. 
The expression for the brightness fluctuations can be fur- 
ther simplified by assuming that the regime T s 3> T ra( j 
holds at the redshift of interest. In this limit, Eq. ([6]) 
reads in Fourier space 

5T h (k,h,e) ~ e" Tc f b (e) [x H A h - x,A Xi - /j'shA,] , 

! (9) 

where /i = k ■ h/k and T b = T^/x^ is the background 
brightness temperature if hydrogen were fully neutral, 
i.e. from Eq. (0} and © 



T h (e) 



%A 10 Ml » b H 2 (l - Y) 
16j4lTO h a 2 ^ t )H(rj t ) 



(10) 



In Eq. ©, we have introduced the perturbed quantity 

kv 

A. u = , 

aH ' 



(11) 



where v is the Fourier transform of the radial baryon ve- 
locity. Assuming that, at EoR, we can neg lect the time 
evolution of the gravitational potentials [39j , one can fur- 
ther approximate A v ~ — A b . 



3. Three-dimensional power spectrum 

The three-dimensional power spectrum of the bright- 
ness temperature fluctuations, PsT h > is defined via 

(5T h (k)ST h (k r )) = (2n) 3 5 3 (k - k')P STb (k) , (12) 

where the brackets denote the ensemble average. One 
can also define 



P(fe) = _P(fe), 



(13) 



such that the isotropic real space variance simplifies to 
(5T£(x)) = [ dlnkV(k), (14) 
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Using Eq. ©, the power spectrum reduces to [44 



P STh (k)=e- 2 ^{xlP h {k) 
+ [x H x i P vi (k) — 



xfP n (k) 



2x m x i P lh {k) 

^xlP v (k)}, 
(15) 

where the (cross) power spectra are defined as in Eq. (fT2|) 
between the perturbation variables Ab, A„ and A Xi . The 
key point of Eq. (1151) is that the /it 4 component depends on 
the baryonic matter power spectrum only since P v ~ Ph, 
which is completely fixed by cosmology. On the contrary, 
the power spectra Pa(k), P v [(k) and Pih(k) are reioniza- 
tion dependent and considered as nuisances. As seen in 
Eq. (lisp , they appear only with an angular dependence 
in or /i 2 , but not in /i 4 . This property is expected 
to play a crucial role in the separation of the cosmologi- 
cal signal from those astrophysical contaminants (75l . l7a | . 
Following these references, we can therefore define 



Po(fc) = xlP h {k) + zfPii - 2x H x i P ib (k), 



P 2 {k) = 2 [x nXl P vi {k) 
P 4 (fc) = x%P v {k), 



c 2 a Pvb{k)] 



(16) 



such that Eq. ([15]) simplifies 



P STb (k) = e- 2 ^T b 2 [P (fc) + fi 2 P 2 (k) + ^Pi(k)] . (17) 

Let us emphasize that the power spectrum Psr b (k) is 
not directly accessible to 21 cm experiments, i.e. one 
cannot directly measure the comoving position r of the 
signal, Fourier dual of the wave vector k. Any experiment 
can however determine angular separations © in the sky, 
and signal frequency differences F around any 21 cm line 
of redshift z. In this paper, we consider a FFTT which 
is designed to map an hemisphere. Following Ref. [iij ]. 
one can divide the sky into small patches for which the 
flat sky approximation is valid. In this limit, and F 
in each patch are directly proportional to the variations 
of the comoving distances components perpendicular r±_ 
and parallel rn to the line of sight with respect to the 
central redshift z. 







D A (z) 



F 



'II 



2/0) 



(18) 



Here Da stands for the angular comoving distance, given 
in a flat Universe by 



(19) 



and y(z) is the conversion factor between comoving dis- 
tances and frequency intervals, 



2/0) 



A 21 (l + z) 2 
H(z) 



(20) 



In Eq. ©, the 21 cm brightness temperature fluctuations 
have been obtained in terms of k while it would be more 
convenient to describe it in terms of Fourier duals of the 



frequency differences and of the angular separations de- 
noted as U|| and u± respectively. Their relation to the 
comoving wave vector k reads 



u± = D A (z)k±, 



'II 



y(z)k\\. 



(21) 



As a consequence, the observable three-dimensional 
power spectrum of the 21 cm brightness temperature is 
measured in it-space and reads 



P S T h (u) 



D 2 A (z)y(z) 



(22) 



4- Angular power spectrum 

Out of the flat sky approximation, one has to derive 
the angular power spectrum of the 21 cm brightness tem- 
perature for each observed redshift z. From an infinitely 
thin shell at redshift z, corresponding to the observed 
21 cm energy e(z), the angular power spectrum is given 
by 



C t {z) 



F e (e,k)Fe(e,k) k 2 dk, 



(23) 



where Fg (e, k) are the multipole components today of the 
21 cm photon distribution function. Their expression 
have been derived in Ref. [42| by solving the Boltzmann 
equation during the dark ages for the perturbed photon 
distribution function (see Appendix [A|l . At EoR, these 
equations remain unchanged but the fractional fluctua- 
tions A H i appearing in Eq. (|A2[) are now given by Eq. . 
The monopole source is still defined as in Eq. (JT)) with a 
background spin temperature that can however be driven 
by the new ionizing sources. Foregrounds could also be 
an extra source of spin temperature fluctuations [66j , but 
we have chosen to consider them as additional nuisance 
power spectra (see Sec. IIIB|) . 

One can recover the usual approximated results by 
keeping only the monopole and redshift distortion terms. 
In that case, the angular power spectrum simplifies to 



Ce(z) ~ Aire 



1 b 



dlnfc \V (k)\ji(&Ve)] 



- V 2 (k)j l (kA Ve )f l '(kA Vt )+V 4 (k)[j'/(kAr h ] 

(24) 

which is the equivalent of Eq. (fT?)) on the full sky. In 
the following, the full sky results are derived by using a 
modified version of the CAMB code [H, [zl, [73| in which 
we have implemented the above-mentioned extra-sources 
associated with the EoR. This allows to integrate the 
complete equations of motion (see Appendix |A"]) under a 
given model of reionization which is discussed in the next 
section. In addition, we consider the same finite width 
Gaussian window functions in frequency as in Ref. (42j, 
such that the observed angular power spectrum is actu- 
ally given by their convolution with Eq. (|23p . Finally, we 
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have cross-checked our computations by comparing the 
modified CAMB results with an independent code directly 
integrating Eq. (f24|) . As expected, both spectra match at 
small angular scales and for small enough window func- 
tion (see Fig. [3]). 



B. Reionization model 

1 . Background 

In order to model the EoR, we consider a toy model 
similar to the one used for CMB data analysis [37| . 
namely a smooth transition from the dark ages ionized 
fraction to complete reionization. A similar model has 
also been considered for 21 cm analysis of the global 
reionization signal in Ref. [13, l7lj . CMB polarization 
measurements allows to infer the value of t, the Comp- 
ton optical depth to reionization defined as 



Vo 



a(r])n e (r])a T dr] . 



(25) 



In this expression, <r T is the Thomson scattering cross 
section while n e stands for the number density of free 
electrons. As explicit in Eqs. (fT7f and (|24[) . Thomson 
scattering at EoR affects the overall amplitude of the 
21 cm power spectra through r c (z), the optical depth to 
the redshift of observation. Let us stress that the latter 
is defined as in Eq. (|25|) but starting at the conformal 
time r] e . Equations (TT71) and ([Mf makes clear that, al- 
ready at the background level, the 21 cm signal at EoR 
is strongly correlated with the total optical depth r that 
can be measured in CMB experiments. In order to make 
the dependency explicit, we have used and extended the 
reionization model used in CAMB [72]. In this model, and 
denoting by x c = n c / (n H i+«.i) the total ionized fraction 2 , 
reionization occurs according to 



1 + /h. 



1 + tanh 



(1 + z rcio ) 3 / 2 - (1 + z) 3 / 2 

Areio 



/ho 

2 



1 + tanh 



A* 



(26) 

where / Ho — 0.08 is the Helium density fraction, assumed 
to become firstly ionized at the same time as hydrogen 3 
and doubly ionized at z — z*. The parameter A re ; fixes 
the duration of reionization while z rc i refers to its central 
redshift. The second term encodes the complete Helium 
ionization and remains always very small. It has been 
kept only for consistency with previous CMB analysis 



2 It can therefore be greater than one when Helium becomes ion- 
ized. 

3 



/He 



m H Y 
m H c 1 - Y ' 



with the same fixed values of = 3.5 and A* = 0.5. 
The last term is a constant encoding the residual ion- 
ized fraction from recombination which does not exceed 
a few thousandth. From Eq. (|2"B1) . one immediately gets 
the associated optical depth by using Eq. (f25j) and con- 
versely. The choice of the power 3/2 in Eq. (|26l) is made 
such that for all A rc i , the corresponding value of r re- 
mains equal to the one that would be associated with an 
instantaneous reionization at z = z ro io (see Ref. [zl|). It 
can be convenient to express A ro ; in terms of redshift 
duration and we will use in the following the parameter 
A z defined by 



A rc 



(27) 



For 21 cm purpose, Eq. (|2"o]l already fixes the evolution 
of the hydrogen neutral and ionized fraction, respectively, 
as 



1 + tanh 



(1 + z rcio ) 3 / 2 - (1 + z) 3 / 2 



(28) 

Concerning the spin temperature, since it is expected to 
be physically driven by ionizing sources, we have chosen 
a similar evolution: 



T s (z) = T 5 



da 



x < 1 + tanh 



(1 + z rcio ) 3 / 2 - (1 + z) 



3/2 



(29) 



The quantity T s da is the spin temperature value during 
the dark ages, just before reionization starts, and it is 
completely fixed by the cosmological model. As a re- 
sult, Eq. (f29| introduces only one extra-parameter, T s max , 
which gives the maximum spin temperature value once 
the universe is completely reionized. Unless specified, we 
assume in the following the fiducial T s max = 10000 K. 

Fig. [3] shows the theoretical three-dimensional and an- 
gular power spectra computed at z = 11 for the back- 
ground cosmological and reionization models just de- 
scribed. The top curve in the bottom panel has been 
obtained by a direct integration of Eq. (|24|) . i.e. including 
only monopole and redshift distortions. As expected, it 
matches in shape and amplitude with three-dimensional 
power spectrum (top panel). The other angular spec- 
tra have been computed with our modified version of the 
CAMB code and are convoluted in frequency with Gaus- 
sian window function of bandwidth B. As discussed in 
Ref. [42j , larger values of the bandwidth damps the over- 
all power but allows the sampling of acoustic oscillations. 
In the limit B — ¥ 0, one recovers Eq. 
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FIG. 3: Three-dimensional power spectra VsT h {k) (top) and 
angular power spectrum CV (bottom) for an observation red- 
shift z = 11 and assuming no foregrounds. The cosmological 
parameters are the ones of Sec. Ill Bl (V = 0.088, z re io — 10.54, 
x H ~ 0.87, T B ~ 1360 K). The C« are convoluted with Gaus- 
sian window functions whose standard deviation is given by 
the bandwidth B. Large band widths damp the overall power 
but render the acoustic oscillations visible [421 ]. The thin shell 
limit has been obtained by integrating Eq. 



2. Foregrounds 

At the perturbative level, the shape and time evolu- 
tion of P;b and P v b and Pi; are expected to depend on 
the reionization details. Following Ref. [44|, we adopt a 
phenomenological parametrization of these power spectra 
as 



Pi(fc) = N a [1 + a Vl fcPn + (fcP u ) 2 ] 7,i/2 P h (k) 
Pb(fc) = N ih exp [-oib fcP ib - (fcP ib ) 2 ] Pb(fc) 



(30) 



where Na, No, are normalization constants and an, aib, 
7ii, Pn and Pib are model parameters to be adjusted to 
match numerical simulations [73, [78[ • Since these spectra 
are nuisance, we moreover approximate P vl {k) ~ — pb(fc) 



TABLE I: Parameters values for the nuisance ionization power 
spectra at z = 11, linearly interpolated from Ref. [3]. 



from Eq. (fTTj). Taking the fiducial cosmological param- 
eters value n h h 2 = 0.023, fl dm h 2 = 0.115, h = 0.71 and 
t = 0.088, the reionization model of the previous sec- 
tion gives z rcio = 10.536 at which the hydrogen neutral 
fraction x u (z re i ) = 0.5. 

In Table HI we give the parameters value associated 
with the reionization power spectra of Refs. 0, [77} at 
2=11. These values will be used in the following to test 
the robustness of the parameter estimation with respect 
to the presence of extra ionizing sources. Although not 
rigorously correct, we assume that the redshift depen- 
dence of the ionizing power spectra is completely given 
by the ionized fractions x m (z) and x t {z), i.e. that the pa- 
rameters of Table |T] remain constant during the beginning 
of the EoR. Notice that the ionized fraction dependency 
is not included in Na and Nr, such that the time evo- 
lution of the ionizing power spectra is well behaved, i.e. 
vanishes as soon as Xi — > 0. 

In Fig. [H we have plotted these spectra for z = 11, 
as well as their effects on the expected three-dimensional 
and angular power spectrum. For this redshift, Xs^iPb 
dominates and reduces the overall power whereas at low 
redshift we would find extra power coming for the domi- 
nation of x 2 Pn. 



C. Fast Fourier Transform Telescope 



In this section we give the ideal characteristics chosen 
for the omniscope used in the following forecasts. The 
FFTT we are considering is in between the original de- 
sign introduced in Ref. @ and a very futuristic concept 
used in Ref. [59j for probing the dark ages. Detecting sig- 
nals from the EoR does not require gigantic array as the 
redshifted 21 cm wavelength from z Te i Q is about a two me- 
ters. However, pushing up the scale of the original FFTT 
concept from 1 km to 10 km would allow to probe the end 
of the dark ages, i.e. up to redshift 20-30. In such a de- 
sign, the sensitivity to dark age signals is quite low and 
forecasts suggest that this would not be competitive with 
CMB (see Sec. IIII B|) . However, including those redshifts 
into a multi-redshifts analysis allows to slice in details 
the transition to the EoR, and as discussed in Sec. IIII CI 
allows a good inference of the background reionization 
properties such as duration or spin temperature evolu- 
tion. 
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no ionizing power 
B =0 

B = 10 MHz 




FIG. 4: Effects of the reionization power spectra at z = 11 
(see Tab. HJ on the three-dimensional power spectrum VsT b 
(top). The bottom panel shows the reionization induced de- 
viations on the angular power spectrum Ce, both for the 
thin shell approximation and a window function of bandwidth 
B = 10 MHz. 



1. Beam, noise and ionosphere 

We consider a squared FFTT design of D = 10 km side 
and constituted of ten million antennas, i.e. one every 
three meters [5J. For simplicity, the beam is approxi- 
mated by a Gaussian whose full width at half maximum 
is fixed to 6f w — 0.89X/D to match the actual square 
aperture. In terms of multipole moments, defining 



4 



7fw 



(31) 



the beam function (normalized to unity for I = 0) reads 

B e = e -^ +1 )/^ . (32) 

For a system antenna temperature given by T sys , even- 
tually including the effects of foregrounds removal, the 



induced flux density fluctuations of one receptor of col- 
lecting area Ai have a root mean square given by ffb\ |8C| 



5S„ = 



271 



sys 



A\ ^t exp B 



(33) 



Here i e xp stands for the exposure time, B is the frequency 
bandwidth and the flux density is related to the bright- 
ness temperature by S v — 2T\,jA\. For aperture synthe- 
sis against iV an t antennas, these fluctuations are reduced 
by a factor ^ N ant (N ant - l)/2 ~ N ant /y/2 correspond- 
ing to the square root of the number of simultaneous 
pair- wise correlations. As a result, the noise fluctuations 
in the brightness temperature for an array of area D 2 
read 



5Tn = * 

V2 N ant 



sys 



(34) 



up to order one factor. The covering factor / cov = 
N ant Ai/D 2 = A cS /D 2 . For a sky map, ST 2 gives the 
expected noise variance in each pixel subtending a solid 
angle (5f2 ~ (A/7)) 2 . The angular noise power is thus 
given by 81] 



sr 2 sn 



X 2 



T 2 



A 2 



T 2 

sys 



2 f2 Uf 
cov JDL cxp 



D 2 l 



^eff /cov Bt, 



(35) 



cxp 



Depending on the scanning strategy, a telescope may 
eventually map a solid angle f2 s ky = 47r/ s k y larger than its 
instantaneous field of view f2f ov . This may be taken into 
account by renormalizing the exposure time by a factor 
fifov/^sky @- The angular noise power finally reads 



C, 



X 2 



T 2 ys 47r/ s 



sky 



sys 



D 2 f ( 



2 Rt 

cov ^"'cxp 



n 



fov 



(36) 



It is well known from CMB data analysis that, for uncor- 
rected signal and noise, the root mean square fluctua- 
tions for the angular power spectrum are given by [82l [83[ 



21+1 



G, 



B 



(37) 



The beam effects can thus be included by defining the 
noise angular power spectrum to be 



(1 

y^f noise w sys 



(38) 



These results also apply to the three-dimensional 
power spectrum. Starting from Eq. (|35[) . the correspond- 
ing noise spectrum p nolsc reads 



P noisc (k) ~ ST 2 SV = C sys D 2 A yB, 



(39) 



where SV = D\y8VLB is the comoving volume of the 
three-dimensional pixel intercepting the solid angle 5VL 
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and for a redshift resolution associated with the band- 
width B [see Eq. (ITS)) ]. From Eq. (|2"2l . one gets the noise 
power spectrum in the u variable 



A 2 T$y$ 47r/ sky 

fcov ^cxp ^fov 



(40) 



Beam effects at small scales can be incorporated by mul- 
tiplying this expression by exp (u 2 j_/£ 2 ). 

Finally, ionospheric absorption becomes important for 
all frequencies approaching the plasma frequency v v ~ 
12 MHz such that atmosphere is opaque to 21 cm signal 
coming from z > 100. For frequencies v 3> v v (z <C 100), 
the ionospheric absorption coefficient Aj ono on the flux 
density can be approximated by [79], 

10 ,o g (A m „)<o. 5 (MM!!;) J . (4i) 

For forecasting, instead of applying this damping directly 
to the expected signal, we prefer to rescale the noise 
power by multiplying Eqs. ([38]) and (j40|) by Af c 



lono 



2. Mock power spectra 

The fiducial design considered in the following has a 
size of D = 10 km, with a covering factor of / cov = 0.1 
and an exposure time of one year, t 



cxp 



3.1 x 10 7 s. 

The system temperature is chosen larger than its typi- 
cal value, T sys = 400 K, to account for foregrounds re- 
moval [12]. For omniscope, the instantaneous field of 
view and sky coverage have been chosen equal to the vis- 
ible sky fi s ky = f2f ov = 27T, i.e. / s k v = 0.5. In Fig. [SJ we 
have plotted the angular power spectrum expected for 
such a design at various observation redshifts z. 



III. FORECASTS ON COSMOLOGICAL AND 
REIONIZATION PARAMETERS 



S io" 



B=10MHz 

rq i | 




B=I MHz 

T| I | 



z=20 
z=30 
z=40 





FIG. 5: Mock angular total power spectra Ce + C° olsc and 
C" olsc in mK 2 for various observation redshifts around reion- 
ization (top) and at the end of the dark ages (bottom). The 
cosmological parameters are the ones of Sec. Ill Bl (r = 0.088, 
z rc i ~ 10.5). The window functions are Gaussian in fre- 
quency space with a respective bandwidth of B = 10 MHz 
(top) and B = 1MHz (bottom). The FFTT design is for 
D = 10 km and allows probing the dark ages for z < 30, above 
which ionospheric absorption damps the signal (see z = 40). 



In order to minimize the parameter space dimensions, 
we consider a simple primordial power spectrum which 
is a power law with unobservable tensor modes. As 
shown below, even for vanishing ionization power spec- 
tra, FFTT alone does not constrain well the optical depth 
t damaging the expected bounds on the background cos- 
mological parameters. As one of the unconstrained pa- 
rameter is the baryon density £lbh 2 , the situation can be 
greatly improved by combining FFTT with CMB data. 
Taking four redshift slices for the FFTT and Planck-like 
typical CMB data, we forecast the expected errors on the 
cosmological parameters and show that all the parame- 
ters entering the reionization models can be inferred. At 
last, the impact of the ionizing power spectra are dis- 
cussed. 



A. Methodology 

In the following, we use both Fisher matrix analysis 
for the three-dimensional power spectra, and MCMC ex- 
ploration of the parameter space for the angular 21 cm 
and CMB power spectra. This allows us to discuss how 
reliable are the expected constraints in presence of de- 
generacies. 

Unless specified, our fiducial cosmological model is 
the one introduced in Sec. [H] complemented with a 
scalar power law primordial power spectrum V s (k) = 
A s (k/K) ns -\ with a pivot scale K = 0.05 Mpc~\ The 
fiducial spectral index has been fixed to n s — 0.973 and 
the amplitude to ln(10 10 A s ) = 3.16. For convenience, we 
have sampled over the CMB parameter 9, the angular 
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scale of the sound horizon at last scattering, instead of 
the Hubble parameter H . This facilitates joint CMB 
and 21 cm analysis discussed later and we have checked 
that this does not affect significantly the forecasts for 21 
cm alone. 



Fisher matrix ana, 



where 



Following the standard approach, given a fiducial set 
of cosmological and reionization parameters {A Q }, the ex- 
pected errors on FFTT measurements of the power spec- 
trum Psr b (u) can be evaluated by using the Fisher ma- 
trix formalism [70]. Let us assume that we can compute 
the likelihood £, i.e. the expected distribution of the 
data given a certain model. The Fisher matrix approach 
assumes that the expected behaviour of the likelihood 
near the maximum is Gaussian such that its curvature 
alone allows to estimate errors on the parameters. The 
fisher information matrix is defined as 



d 2 ln£ 

dX a dX b 



(42) 



where the brackets denote an ensemble average. When 
all the parameters are estimated simultaneously, the 
marginalized error on a given parameter reads a\ a > 



F ) aa - The Cramer-Rao inequality (">") empha- 
sizes the fact that the Fisher matrix approach always 
gives a local optimistic estimate of the errors. It becomes 
an equality when the likelihood is a Gaussian around its 
maximum. In the following, we estimate the expected 
variance on the measured parameters A a from the diag- 
onal elements of the covariance matrix, i.e. from the 
inverse Fisher matrix elements 



(43) 



For the purpose of the Fisher matrix analysis, as al- 
ready mentioned in Sec. Ill A 31 one can divide the sky 
into small patches for which the flat sky approximation 
is valid. Given an angular patch O and a frequency bin 
of size F, the probed comoving volume is 



V 



2 D 2 A Fy. 



(44) 



Following [44(, one can consider step sizes in the dual 
it-space, as follows 



2?r 2tt 

-L = -q i Su \\ = -y 



(45) 



where 8 is taken to be lower than typically 1 rad and F 
is set by the frequency size of the redshift bin. Assuming 
that PsT h {u) is Gaussian distributed, we can approxi- 
mate the Fisher matrix by 



2 ^ 

Mil ,111 



Nr. 



dPsrJu) dP 5Tb (u) 



[P STb (u) + P» 



d\ a 



(46) 



(47) 



is the number of independent cells probed 4 for a given 
value of fc, or equivalently of it. Notice that given 
Eqs. (gl|) and (|35]1. iV c reduces to 



Nr. 



47r/ sk y 2nk± 
Q 2 6k± 



(48) 



In the following, we will be interested in comparing the 
results of the Fisher analysis on the three-dimensional 
power spectrum PsT b (u) to the results of a MCMC anal- 
ysis which makes use of the angular power spectrum Q. 
From our experimental design, the measured 21 cm power 
spectrum is however convoluted in frequency by a Gaus- 
sian window function of width B. As the result, and un- 
less specified, we will also compute the Fisher matrices of 
Eq. (I46[) from the three dimensional power spectrum at 
feii = after a convolution with the same Gaussian win- 
dow function. Compared to Ref. (44|, this allows us to 
consider larger redshift bins, and thus larger bandwidths, 
in which the power spectrum and background quanti- 
ties are allowed to vary. In Ref. [44j, the FFTT was 
assumed to have a small bandwidth of B ~ 0.15 MHz. 
Around each observation redshift, the signal on a bin of 
size 6z ~ 0.5 is expanded in Fourier modes u\\, up to 
«™ = 2-k/B. Within our approach, this Fourier expan- 
sion is now replaced by a convolution of the signal with 
the frequency window function associated with the mean 
redshift of observation. Combining multiple redshifts is 
still done by adding up their respective Fisher matrices. 

Let us stress that when the bandwidth is small enough, 
there is no need for the convolution. Summing the Fisher 
on u\\ running on positive and negative values with |u™ ax | 
set by the frequency resolution of the experiment, i.e. the 
bandwidth B, gives the same result than subdividing the 
redshift bin in an identical number of slices to the number 
of parallel modes. We have cross-checked our method by 
using the FFTT specifications and reionization model of 
44( | . for which we recover the same results, 



Ref. 



the 



small bandwidth limit. 

Concerning the transverse modes, since the volume of 
the experiment is finite, only a discrete number of modes 
is theoretically accessible to the analysis. In Eq. (|46p . the 
sum over u± runs over positive values from u™ 111 = 27r/0, 
to w™ ax ~ 2nD/X, the angular resolution being set by the 
longest baseline of the experiment, D. 

Finally, let us also mention that when constructing the 
Fisher matrix, one has to pay a special attention to the 
numerical accuracy at which the power spectrum deriva- 
tives have been estimated. Indeed, the evolution of the 



4 The factor 1/2 in Eq. (146 I t takes into account the fact that k and 
— k are not independent. 
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21 cm power spectrum with the cosmological parame- 
ters can be highly non linear and is furthermore red- 
shift dependant. As discussed in Ref. [84|, one has to 
find a good compromise between small step sizes in the 
numerical evaluation of the derivatives and the appear- 
ance of numerical noise. For each parameter, a step size 
corresponding roughly to 2% of the fiducial values was 
found to be a good compromise. For each redshift, the 
power spectrum needs to be calculated for a large num- 
ber of wavelength modes (typically 3000) and, in order to 
spare computing time, it has been interpolated using 20 
/c-values per unit logarithmic interval. The derivatives 
have been evaluated by means of a four-point method 
for well behaved parameters whereas for n s , A s , plus the 
nuisance reionization parameters Na, Rn, aii,7u, Mb, -Rib 
and ait,, the derivatives have been determined analyti- 
cally. Finally, to keep the numerical noise under control, 
the transfer functions have been computed using high 
precision settings within CAMB. 



2. Monte-Carlo-Markov-Chains 

From the mock data associated with the FFTT design 
of Sec. Ill Cl forecasts can be derived through MCMC ex- 
ploration of the parameter space provided one specifies 
the likelihood. As for the Fisher analysis, the mock data 
{Ci} are assumed to be associated with a set of fidu- 
cial parameters {A a }. For a full sky analysis, one can 
show that the sampling distribution followed by the Ce 



is a gamma distribution 



For a cut sky, with an 



isotropic beam and noise, one can use the approximated 
likelihood EE H3 



Here Di denotes the theoretical angular power spectrum, 
plus noise, i.e. 



D. = Ce 



Cf 



B 



2 ' 



(50) 



associated with any tested set of parameters {A a }. Simi- 
larly, Di are the equivalent for the fiducial Cg. From our 
modified version of the CAMB code, and FFTT specifica- 
tions, the Di can be computed for any input value of the 
cosmological and reionization parameters {A a }. Using 
MCMC sampling with the likelihood of Eq. @ allows 
to extract their posterior probability distribution for a 
given fiducial model {A a }. In the following, we have used 
a modified version of the public code CosmoMC coupled 
with our modified CAMB code l73Tl. 
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FIG. 6: Expected marginalized posterior probability distri- 
butions (solid) and mean likelihood (dotted) for 10 km FFTT 
data alone, at z = 30 (dark ages). The total optical depth r, 
primordial amplitude As, baryon density Qhh 2 are completely 
degenerated and not constrained (see text and Fig. . 



B. Single redshift from dark ages 

In this section, the above-described MCMC analysis 
is applied to one dark ages redshift z = 30. As dis- 
cussed in Sec. Ill Cl the 21 cm physics during dark ages is 
free of reionization uncertainties, up to the optical depth. 
Within our design, it has a low signal-to-noise ratio whose 
only interest is to exacerbate the degeneracies induced by 
the optical depth into the other cosmological parameters. 
In Fig. [SI we have plotted the marginalized probability 
distribution on the cosmological parameters for the same 
reionization model used for CMB data analysis, namely 
only t (or z reio ) is varying. 

Both the primordial power spectrum amplitude A s and 
r enters into the amplitude of the 21 cm signal and one 
expects them to be degenerated. As Fig. [5] and Fig. [7] 
show, these parameters are in fact mostly degenerated 
with Q\>h 2 . This is expected because, as can be seen in 
Eq. (I10[) . the baryon density also enters into the over- 
all amplitude of the signal. Moreover, even for a fixed 
r, and fixed observation redshift z, r c still inherits ex- 
tra dependencies on the other cosmological parameters 
from their influence on the background evolution. As a 
result, all parameters end up being more or less contam- 
inated by the reionization parameter r, which seriously 
damages the expected constraints compared to a situa- 
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FIG. 7: One and two sigma confidence interval of the two- 
dimensional marginalized posterior probability distributions 
(solid lines) for a 10 km FFTT design at z = 30. The shading 
traces the two-dimensional mean likelihood. 



tion in which r would be perfectly known. The notable 
exception concerns the spectral index n s , simply because 
the number of modes accessible to the FFTT is large 
enough to detect its effect over background degeneracies, 
these ones affecting all modes in a smoother way. Nev- 
ertheless, the forecasts on n s do not compete with CMB 
bounds. Let us stress that we have used only one redshift 
slice to emphasize the above-mentioned degeneracies and 
one could add other redshifts to improve the constraints. 
However, for the dark ages, the signal dominates over 
the noise only around z ~ 30 (see Fig. and the above 
degeneracies still remain. The situation is better at EoR. 



C. Epoch of reionization 

The fiducial model is unchanged but the observation 
redshifts now probe the beginning of the EoR. 



1. Optical depth induced degeneracies for one redshift 

As for the dark ages, we have first performed our analy- 
sis for a single redshift slice z = 1 1 , close to the maximum 
signal-to-noise ratio of the 10 km FFTT and sampled at 
B = 10 MHz. In an optimal situation, the ionizing power 
spectrum can be switched off and the background reion- 
ization model has only r varying. 

The MCMC forecasts are represented in Fig. [8J and 
Fig.|H]as one and two-dimensional posterior distributions. 

As for the dark ages, the degeneracies induced by r 
are still present and propagate to the other parameters 
as well. However, r, A s and fl^h 2 are now constrained. 
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FIG. 8: Expected marginalized posterior probability distri- 
butions (solid) and mean likelihood (dotted) from a single 
redshift observation at z = 11 (reionization). The non- 
Gaussian posteriors and "shoulders" apparent for some pa- 
rameters comes from the large correlations induced by r (see 
Fig. [5}. 



This comes from the signal to noise ratio at z = 1 1 which 
is much larger than for the dark ages, but also from the 
stronger effects induced by r which makes the degenera- 
cies non-linear. As seen in Fig. [51 the marginalized prob- 
ability are non- Gaussian and distorted by their correla- 
tions. This is particularly visible in the two-dimensional 
likelihood in the plane (r, Ob/i 2 ) (see Fig. EJ) - Compared 
to the dark ages, varying r affects the background sig- 
nal amplitude at EoR through x a and T s . Compared to 
previous works derived with Fisher matrix methods, our 
forecasts are slightly less optimistic, but as it should be 
clear from Fig. [9j not including r reduces the expected 
error bars on various other cosmological parameters (see 
below). The need of breaking degeneracies clearly sug- 
gests the inclusion of additional data, either 21 cm from 
other redshift slices, or CMB data. 

We have also performed a Fisher matrix analysis under 
the same model and experimental hypothesis, both on the 
convoluted and non-convoluted 21 cm three-dimensional 
power spectra. Our results are reported in Tab.HIl With- 
out convolution, the strong degeneracy between r and 
\nA s leads to a ill-conditioned Fisher matrix that can 
be hardly inverted. This can therefore induce important 
numerical errors in the forecasts. Physically, this means 
that one combination of parameters is completely un- 
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FIG. 9: Two-dimensional posteriors for a single observation 
redshift at z = 11. The solid lines represent the one and two- 
sigma confidence interval while the shading traces the two- 
dimensional mean likelihood. Notice the strong degeneracies 
induced by r which boost the expected variance on Qbh 2 for 
instance. 



constrained from the Fisher point of view: the Gaussian 
approximation would predict some flat directions for the 
likelihood. In such a situation, Fisher results can still be 
derived provided the problematic parameters are artifi- 
cially removed to recover a well-behaved matrix. 

Considering the convolution in frequency reduces the 
degeneracy between the parameters r and In A s but still 
does not cure the existence of flat directions. In Tab. [Ill 
we have reported the forecasts when removing either 
ln^4 s or t from the Fisher analysis. Forecast on the 
other parameters can be slightly over-pessimistic com- 
pared to the MCMC. This is expected because, even 
with the convolution included, the Fisher results still ap- 
proximate any degeneracies as linear. Fig. |H] shows that 
this is clearly not the for instance in the plane 

(flbh 2 , r) where the two-sigma contour is actually multi- 
valued. Let us notice that removing r in the Fisher ma- 
trix yields an order of magnitude wrong prediction for 
In A s , whereas removing In A s gives a standard deviation 
for r compatible with the MCMC result (see Tab. |TT]). 
These observations reinforce the need of MCMC meth- 
ods as soon as some degeneracies are not under control. 
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TABLE II: Standard deviations on cosmological parameters 
for FFTT-only single redshift measurements at z = 11 from 
both MCMC method (M) and Fisher matrix analysis (F). 
The Fisher results are given for both the non-convoluted and 
convoluted power spectrum and by removing from the anal- 
ysis one of the strongly degenerated parameters, either r or 
ln(10 10 A s ). 



2. Multi-redshifts tomography 

Still with vanishing ionizing power spectra and a reion- 
ization model having only r varying, we can reduce the 
FFTT degeneracies by adding extra redshift slices. As 
before, we have performed MCMC analysis on the full 
sky mock data and Fisher matrix methods on the asso- 
ciated 21 cm power spectrum. 

In Fig. [I0]and Fig.QUwe have represented the MCMC 
posteriors obtained by considering two redshift slices, z = 
11, z — 12, and four redshift slices z = 10, z = 11, z = 12, 
z = 13. The telescope is, as before, the 10 km FFTT with 
B = 10 MHz and all other parameters are the same. 

Compared to z = 11 alone, the degeneracies have been 
reduced. In fact, one can check that a data analysis per- 
formed at z = 12 alone would produce similar posteri- 
ors as those obtained for z — 11 with however slightly 
different degeneracy directions in the parameter space. 
Combining both can be viewed as keeping only the in- 
tercepting probability contours of each redshift slice. For 
r, the standard deviation ends up being reduced by al- 
most a factor of two, while this is not as significant for 
the other parameters. Their posteriors do no longer ex- 
hibit "shoulders" or evident non-Gaussian shapes, how- 
ever Fig. [TT] shows that they still are correlated along 
some peculiar directions. Adding more redshifts again 
improves the situation, but the correlation trend does not 
disappear due to the common physical origin of all 21 cm 
signals at reionization. The inner red curves in Fig llOl and 
Fig [11] are the same analysis performed on four redshift 
slices, scanning the whole reionization duration: z = 10, 
z = 11, z = 12 and z = 13. The standard deviations 
are slightly reduced, but less than the expected \/2 fac- 
tor because both z = 10 and z — 13 correspond to a low 
emission signal (see Fig. [S])- Let us notice that a band- 



14 









0.021 0.022 0.023 0.024 0.025 


0.112 0.114 0.116 0.118 








.038 1.04 1.042 1.044 1.046 1.048 

e 




0.97 0.972 0.974 0.976 

n 








0.084 0.086 0.088 0.09 0.092 0.094 

X 


3.05 3.1 3.15 3.2 3.25 3.3 
ln[10 10 AJ 









70 71 72 73 0.265 0.27 0.275 0.28 

H Q » 



FIG. 10: Marginalized posteriors (solid) and mean likelihood 
(dotted) from two and four observation redshifts at reioniza- 
tion. The black outer curves are for two redshift slices z = 11, 
z = 12 while the red inner ones for four redshift z — 10, 
z = 11, z = 12 and z = 13 (see also Fig. [TT|) . 



width B = 10 MHz corresponds in term of redshift to a 
width Sz ~ 1 (at z = 11). Requiring more redshift slices 
at reionization would therefore also demands a smaller 
bandwidth. 

Let us turn now to the Fisher analysis. The forecasts 
for the four redshifts z = 10, z = 11, z = 12 and z = 13 
and the convoluted 21 cm power spectrum are given in 
Tab. IIII1 Our results agree relatively well with those from 
the MCMC method, even if the forecasts are found to be 
about 20% more stringent with the Fisher matrix anal- 
ysis. We have also checked that the degeneracy direc- 
tions between parameters were matching the ones ob- 
tained from the MCMC and represented in Fig. [TT1 



3. Adding Planck-like CMB data 

The previous results call for the incorporation of other 
cosmological data having different correlations. In the 
following, we incorporate CMB data and consider a 
typical Planck-like experiment. For our purpose, we 
have generated mock CMB data exactly as described in 
Sec. Ill CI but for a full sky CMB experiment (/ s k y = 1), 
with a Gaussian beam of resolution 9{ w = 7' and a con- 
stant noise for each pixel given by (ST?)= 2 x 10~ 4 /iK 2 
(twice for the E'-mode polarization) [88], |89[ . Perform- 




0.021 , 0.024 0.021 , 0.024 
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FIG. 11: Two-dimensional posteriors obtained by combining 
the FFTT data from multiple redshift slices. The shading 
traces the mean likelihood associated with the combination of 
two redshifts z = 11 and z = 12. The black largest contours 
are the associated one- and two-sigma confidence intervals 
whereas the smaller red contours are those associated with 
combined data for four redshifts (z = 10, z = 11, z — 12 
and z = 13). The degeneracies induced by r are accordingly 
reduced. 



ing a MCMC analysis with Planck-like mock data alone, 
and for the same fiducial parameter values, yields the 
posterior plotted in Fig. Q21 To ease the comparison, we 
have reported on the same figure the two redshifts FFTT 
forecasts. 

CMB data give much stronger constraints on ft^h 2 , 
A s and 9 than the considered 10 km FFTT design. This 
is expected because the overall amplitude of the CMB 
signal depends mainly on A s and not fl\>h 2 contrary to 
the 21 cm brightness. Concerning 0, it is an optimal pa- 
rameter for CMB by design whereas Hq shows similar 
posteriors for both CMB and 21 cm. Interestingly, the 
spectral index is, by a factor of two, more constrained by 
21 cm. One can also check in Fig. [3] that this would also 
the case with only one observation redshift at z = 11. 
This is again due to the huge level arm provided by the 
FFTT resolution on the power spectrum shape. More- 
over, Fig. [T2l shows that, in spite of the degeneracies, the 
marginalized posterior for r remains sharper for FFTT 
than for CMB. This could wrongly suggest that adding 
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FIG. 12: Typical marginalized posterior distributions for a 
Planck-like CMB experiment (black solid lines). The dashed 
blue curves are the two redshift forecasts for FFTT (same as 
in Fig. HOl) . 



CMB data may not help to improve the FFTT con- 
straints over r. However, r is mostly correlated with 
Q,\,h 2 which is well bounded by CMB alone. As a result, 
CMB data are expected to indirectly improves the FFTT 
accuracy on r. 

We have checked that this is indeed the case by per- 
forming a MCMC analysis based on Planck-like CMB 
data together with the four redshifts FFTT data. The re- 
sulting posterior probability distributions are represented 
in Fig Q21 The expected standard deviation on r is re- 
duced by an order of magnitude. In Fig. [TU the two- 
dimensional posterior confirms that such an improvement 
comes from the significant reduction in the variance of 
il^h 2 due to CMB data. For comparison, we have also 
reported in this figure the one- and two-sigma confidence 
intervals for the four redshifts FFTT alone (see Fig. ITT1) . 
For some parameters, the improvement is such that they 
may lie outside the displayed plane. 

The forecasts derived from the Fisher matrix are re- 
ported in Tabic IIIII As the degeneracies are even more 
reduced by the inclusion of CMB data, the Fisher re- 
sults agree well with those from the MCMC analysis. 
The degeneracy axes associated with one-sigma ellipses 
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FIG. 13: Marginalized posteriors from the combination of the 
FFTT four redshifts data (z = 10, z = 11, z = 12 and z = 13) 
with Planck-like CMB data. The reduced variance on fl^h 2 
coming from CMB data kills most of the degeneracies asso- 
ciated with the 21 cm signal. Four redshifts FFTT without 
CMB is represented in Fig. [TU] 



are also found to be in agreement with the MCMC two- 
dimensional posteriors (not represented). 

As the Fisher matrix approach can be trusted in this 
regime, one can discuss the effect of increasing the num- 
ber of modes by reducing the bandwidth. For B = 
0.1 MHz, still using the same four mean redshifts, with a 
bin size of 5z = 0.5, it is no longer necessary to include 
the convolution. In that case, there are about 40 parallel 
modes accessible along the line of sight, up to a non-linear 
cut-off fixed at k\\ = 2Mpc _1 0] (see Tab.|m|. 

In the following sections, we relax some of the reioniza- 
tion assumptions and consider that duration and maxi- 
mal spin temperature are no longer known. 

4- Spin temperature evolution and reionization duration 

The reionization model of the previous sections was 
involving an unique varying parameter r. One may now 
wonder how the forecasts are affected by the reionization 
duration A z and the maximal spin temperature T s max . 
Again, the fiducial model in unchanged, i.e. has A 2 = 0.5 
and T™ ax = 10000 K. 

Fig. Qj)] and Fig. [If)] show the probability distributions 
obtained from the MCMC exploration based on the same 
data, i.e. four redshifts FFTT plus CMB Planck-like 
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Method 


fl h h 2 




8 


n s 


T 


ln(10 10 A s ) 


Ho 


FFTT alone 


M 
F 

F (B = 0.1MHz, no conv.) 


5.4 x 10" 4 

4.2 x 10" 4 

2.3 x 10" 5 


1.2 x 10" 3 

9.5 x 10" 4 

1.6 x 10" 4 


1.6 x 10" 3 

1.3 x 10" 3 

1.4 x 10" 4 


1.3 x 10" 3 
1.0 x 10" 3 
1.7 x 10" 4 


1.5 x 10" 3 
1.2 x 10" 3 
4.9 x 10" 5 


3.8 x 10" 2 

2.9 x 10" 2 
2.3 x 10" 3 


0.55 
0.45 
4.9 x 10" 2 


FFTT + CMB 


M 
F 

F (B = 0.1MHz, no conv.) 


3.0 x 10" 5 
2.8 x 10 -5 

8.1 x 10" 6 


2.9 x 10" 4 
2.6 x 10 -4 
3.9 x 10" 5 


2.0 x 10" 4 

2.2 x 10" 4 

8.3 x 10" 5 


8.6 x 10 -4 
7.2 x 10" 4 

5.7 x 10" 5 


1.8 x 10" 4 
1.6 x 10 -4 
3.6 x 10" 5 


9.5 x 10 -4 
9.0 x 10" 4 
3.4 x 10" 4 


0.15 
0.13 
1.2 x 10" 2 



TABLE III: Expected standard deviations for the basic parameters in the case of four redshifts FFTT alone and complemented 
with Planck-like CMB data. The reionization model has only one varying parameter r, the ionizing power spectra are assumed 
to be negligible (for the MCMC method, Ho is a derived parameter and is reported in the last column). The method refers to 
either MCMC (M) full sky analysis or the Fisher matrix (F) on the 21cm three-dimensional power spectrum (with convolution). 
The third and sixth lines show the Fisher results (without convolution) in the small bandwidth limit B = 0.1 MHz, with Sz = 0.5, 
for which 40 wy modes would be accessible. 
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FIG. 14: Two-dimensional marginalized probability distribu- 
tion for the combined four redshifts FFTT and Planck-like 
CMB data. The inner black solid lines are the one- and 
two-sigma confidence intervals, the shading traces the mean 
likelihood. The outer red lines are the one- and two-sigma 
confidence intervals for four redshifts FFTT data alone. The 
precise determination of fib ft 2 from CMB strongly reduces the 
(r, Q.bh 2 ) correlations associated with the 21 cm signal. 



data. The two additional reionization parameters are 
very well constrained. Even for the large chosen fidu- 
cial value T s max , the expected standard deviation do not 
exceed 20 K. This directly comes from the FFTT tomog- 
raphy. Along the four redshifts probed, the 21cm signal 
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FIG. 15: Posterior probability distributions from combined 
FFTT four redshifts and Planck-like CMB data. The back- 
ground reionization model has three varying parameters, r, 
A z and T, max , all well reconstructed. 
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TABLE IV: Standard deviations for the combined four red- 
shifts FFTT and CMB for a reionization model with three 
varying parameters r, A z and T s max . There is no loss of accu- 
racy compared to Table. [Til] showing that redshift tomography 
allows the complete reconstruction of the reionization param- 
eters. The label "M" and "F" respectively refers to MCMC 
and Fisher matrix analysis (with convolution). The Fisher 
matrix has to be regularized at the expense of loosing infor- 
mation by removing one of the degenerated parameters, here 

jimax 



pear for A z . This is again not surprising because those 
parameters are correlated (see Fig. ITTJl) . 

In the next section, we relax the assumption that the 
ionizing power spectra are negligible and consider their 
impact on the model parameter forecasts. 



FIG. 16: Two-dimensional posterior probability distributions 
from combined FFTT four redshifts and Planck-like CMB 
data (see also Fig. I15[) . The background reionization model 
parameters, t, A z and T s max , are weakly correlated and well 
reconstructed. 



is very sensitive to the neutral hydrogen fraction x H (z) 
and the value of T a {z) thereby providing accurate fore- 
casts for the related parameters. As visible in Fig. [T51 
ymax - g gijg^iy correlated with A z but this is expected 
in view of our smooth reionization model (see Fig.Q]). In- 
creasing T s max at a given redshift mimics a reduction in 
A z . The correlations remain under control because we 
use more than one redshift. These results also suggest 
that taking the asymptotic limit T s — > oo at reionization 
is not always justified and throws away some potentially 
available information. 

We find that the Fisher matrix analysis, with convolu- 
tion included, does not resolve completely the degeneracy 
between T s max and A 2 thereby rendering the Fisher ma- 
trix singular. As before, one needs to remove the incrim- 
inated flat directions, and we had to pull out T" lax from 
the analysis 5 . In this case, the predicted standard devi- 
ations are in good agreement with the MCMC results, 
whereas larger differences between the two methods ap- 



5 Removing A z still yields an almost singular matrix generating 
spurious numerical errors at inversion. 



5. Nuisance ionization power spectra 

As in Ref. [lij], we consider a less ideal situation in 
which the ionizing sources affect the 21 cm perturba- 
tions by switching on the extra power spectra P;b, P& 
of Eq. ([TT)|) . According to the discussion of Sec. Ill B 21 
the redshift evolution of our toy ionizing spectra is com- 
pletely encoded into x H (z) and x^z); as in Eq. (fT6| . In 
total, this adds seven nuisance parameters to the analy- 
sis, namely Na, R lu an, jn, N^,, and a^- In a more 
realistic situation, the ionizing spectra may involve more 
parameters and have certainly more complicated shapes, 
but their redshift evolution should still be correlated with 
the background reionization history. 

Using again the four- redshifts FFTT plus Planck-likc 
CMB data, we have plotted in Fig. [17] and [18] the re- 
sulting marginalized probability distributions when the 
ionizing spectra are switched on. The fiducial values of 
all parameters are the same as before while the seven nui- 
sance parameters have been fixed to the fiducial values 
of Table H 

The most important result is that the posteriors for 
the cosmological parameters are almost unchanged com- 
pared to Fig. [TS] thereby confirming that the presence of 
nuisance ionizing spectra should not prevent cosmolog- 
ical parameter estimation [Hj]. The expected standard 
deviations have been reported in Table. [V] and, up to 
an increase for Vt\yh 2 and r, the bounds of all the other 
cosmological parameters are identical to the ones of Ta- 
ble. IIVI Let us stress again that this result incorporates 
marginalization over all background reionization histories 
accessible within our model. 
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FIG. 17: Posterior probability distributions in presence of the ionizing power spectra. The background reionization model has 
three varying parameters, r, A z and T s max while the ionizing spectra involves seven nuisance parameters Nu, an, jn, Rn, ATjb, 
«ib and Rn,. The posteriors of Fig. [15] have been reported as dashed curves to ease comparison. 



The situation is different for the background reioniza- 
tion parameters. The standard deviation associated with 
the posteriors of A z and T s max is more than one order of 
magnitude larger than those obtained without ionizing 
spectra. As can be seen in Fig. [18j the two-dimensional 
posteriors exhibit strong correlations between the nui- 
sance parameters and A z and T s max . This can be un- 
derstood from Fig. [4j Around the redshift of maximum 
sensitivity (z ~ 11), the total ionizing spectra is dom- 
inated by Pit,. Its amplitude N-^ therefore affects the 
overall observable signal and should therefore correlate 
to all the other amplitude related parameters. However, 
since we use CMB data, the parameters flhh 2 and r are 
relatively well constrained, independently of the values 
of iVit,. They are therefore weakly affected by these new 
correlations. It remains the neutral fraction x H , which 
is precisely given by the background reionization his- 
tory, i.e. by the value of A z and T s max . One can see in 



Fig. [18] that the two-dimensional posteriors in the planes 
(A z ,Nib) and (T s max , iVjb) are indeed degenerated and 
even multi-valued. The "strips" of maximum posterior 
values are due to the presence of the acoustic oscillations 
which favour only some peculiar combination of ionizing 
amplitude and neutral fraction. 

When marginalizing over the cosmological and back- 
ground reionization parameters, Fig.[T71and Tab. IVl show 
that the nuisance parameters iVjb, c^b and Rib can actu- 
ally be inferred. This is not surprising as they encode 
the dominant contribution for the ionizing power spec- 
tra. On the other hand, although the posteriors of Nn, 
Ru, an and jn peaks at some preferred values, these are 
significantly biased compared to the fiducial ones. Again, 
this can be understood from Fig. [4] The Pa contribution 
only shows up at small scales, typically £ > 10000, where 
the noise power spectrum becomes of comparable ampli- 
tude to the signal (see Fig. [3]). 
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FIG. 18: Two-dimensional posterior probability distributions 
in presence of the ionizing power spectra. The correlations 
induced by the nuisance parameter Nih severely damages the 
expected constraints on the background reionization parame- 
ters A z and T s max . The multi- valued confidence intervals are 
the results of the acoustic oscillations (see text). 



Unsurprisingly, the above mentioned degeneracies are 
problematic for a Fisher matrix analysis. As before, the 
only solution to render the Fisher matrix regular was to 
remove the parameter T s max , whereas removing A z would 
not fix the ill-conditioning. As can be seen in Tab.fVl the 
Fisher method (convolution included) has difficulty to re- 
produce all of the MCMC results accurately. For some 
parameters, precisely the ones which remains weakly sen- 
sitive to the ionization parameters, the Fisher forecasts 
match well with the MCMC. This is the case for the 
standard cosmological parameters (up to il^h 2 ). For the 
background reionization parameters, as for instance A z , 
the Fisher expected variances can be out by a factor five, 
which slightly bias the expected variance of r and fl^h 2 . 
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TABLE V: Standard deviations in presence of the ionizing 
power spectra. The mock data are, as before, four redshifts 
FFTT and Planck-like CMB. Although there is no loss of 
accuracy compared to Table. ITVl for the cosmological parame- 
ters, the expected constraints on the background reionization 
parameter A z and T" s max are severely damaged due to the ap- 
pearance of new degeneracies (see Fig. [THl) . 



Turning off T s max also kills the multi-modal behaviour 
of the likelihood and the Fisher results end up luckily 
giving the correct results for the ionizing power spectra 
amplitudes, Na and iVib- 

To summarize, we have found that even in presence of 
ionizing power spectra, combining multi-redshift FFTT 
data with CMB data can still be used to reconstruct 
the background reionization history together with the 
more usual set of cosmological parameters. The expected 
variances of the reionization parameters are nevertheless 
severely damaged compared to no ionizing sources. 



IV. CONCLUSION 

In this paper, we have quantitatively shown that om- 
niscopes could be used to constrain the background reion- 
ization history while being not sensitive to the zero mode 
of the brightness temperature. 

For this purpose, we have considered a simple, but con- 
sistent, reionization model completely determined by the 
total optical depth t, the reionization duration A z and 
the asymptotic spin temperature T s max . We have used 
both Fisher matrix approaches on the three-dimensional 
power spectrum and MCMC methods on full sky simu- 
lated data to forecast the expected variance of all model 
parameters. Our results suggest that it is crucial to com- 
bine multiple redshifts at the EoR with CMB data to 
keep all degeneracies under control, in particular owing 
to an accurate determination of ft^h 2 . When this is not 
the case, we have shown that the Fisher predictions can 
be quite inaccurate, eventually being over-optimistic, but 
also over-pessimistic when the conditions for the Cramer- 
Rao inequalities are no longer met. For those situations, 
only the MCMC methods ends up being usable. 

Within combined data, the perspectives are quite good 
as the background history can be fully reconstructed, 
even in the presence of correlated nuisance ionizing 
sources. Although we have discussed only one kind of 
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FFTT telescopes, our results could be easily generalized 
to other configurations by a proper rescaling of the res- 
olution and noise parameters of Sec. HI1 Still, our work 
would need to be extended with a more complete reion- 
ization model, eventually adjusted to numerical simula- 
tions in order to include a realistic dependence between 
the reionization parameters and the redshift evolution of 
the ionization power spectra. Another extension would 
be the inclusion of exotic reionization sources, such as a 
decaying or annihilating dark matter component. 
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Appendix A: Boltzmann equations 

In this appendix, we recap some of the equations driv- 
ing the dynamics of the 21 cm brightness fluctuations as 
they are derived in Ref. (42J. Denoting by f(rj, e) the 
photon distribution function of 21 cm photons, where e 
is the redshifted energy of 21 cm photons (7b = e//2 in 
the Rayleigh- Jeans approximation), the multipole com- 
ponents of the photon distribution function Fi(r/, e, k) are 
defined by 

/d fr 3 
— ^ (-i) l F t Yr(k)Y t m (*), 

(Al) 

where the Y™ are the spherical harmonics. From the per- 
turbed Boltzmann equation, they are found to verify [42| 
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In this equation, the conformal distance to the redshift of 
observation is given byA?7 e = rj e — 770, r c is the Thomson 
scattering optical depth, r e the optical depth of 21 cm 
photons and r e = r e e _rt /(l — e~ Te ) encodes the amount 
of rescattered 21 cm photons. The Legendre polynomi- 
als have been denoted Pi while Of stand for the multi- 
pole moments of the CMB temperature anisotropies [13] ■ 
Both f e and g e are function of e only conveniently defined 
from the background distribution function 



1 _ g-T2l(»?) 

1 - e~ T > 



(A3) 



and its derivative 



e5 e /(r?,e) = g, 



l-e- T "W / e (d v T 21 )t 



-T 2 l(v) 



1-e- 



U(r] e ) 1 - e— ■ 



In the above equations, 721(77) = T e H(?7 — r) e ), H(x) being 
the Heaviside step function. As discussed in Sec. Ill A 41 
the the small angle limit, assuming r c -C 1 such that r e ~ 
1, neglecting photon perturbations, one recover Eq. © 
but for 5f and the multipole moments read 



F e ~ e Tc T h [(x H A h - XiA^) j e (kAr] e ) 
+ x n A v j' e '(kAr] e )] . 



(A5) 



(A4) Plugging this equation into Eq. $23j) gives Eq. (|24|) . 
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